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We predict a new mechanism of enhancement of ferromagnetic phase transition temperature Tc 
in uniaxially stressed diluted magnetic semiconductors (DMS) of p-type. Our prediction is based on 
comparative studies of both Heisenberg (inherent to undistorted DMS with cubic lattice) and Ising 
(which can be applied to strongly enough stressed DMS) models in a random field approximation 
permitting to take into account the spatial inhomogeneity of spin-spin interaction. Our calculations 
of phase diagrams show that area of parameters for existence of DMS-ferromagnetism in Ising model 
is much larger than that in Heisenberg model. 
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I. INTRODUCTION 

Recently (see^-^ and references therein), great ad- 
vances have been made in the problem of ferromag- 
netism (FMm) of p-doped diluted magnetic semiconduc- 
tors (DMS). The question about influence of different 
physical phenomena in DMS on the critical temperature 
Tc of ferromagnetic (FM) phase transition is of prime in- 
terest for searching for the future study trends. After 
Pashitskii and Ryabchenko prediction of FMm in DMS,'^ 
the competition between FM correlations mediated by 
indirect long-range spin-spin interaction and direct short- 
range antiferromagnetic (AFM) interaction is considered 
to be decisive in the formation of FM state in DMS. In 
other words, to obtain FM state with high enough Tc, we 
need to inhibit the AFM contribution. 

Due to short-range character of AFM interaction, only 
close pairs of magnetic ions contribute to it. This con- 
tribution can be minimized by decreasing the magnetic 
ion concentration n^. On the other hand, when rii is too 
small, ferromagnetism can be destroyed. To retain FM 
ordering in this case, we may increase the carriers con- 
centration Uc (in p-type DMS this corresponds to holes 
concentration ). The calculations of Tc performed in 
a mean field approximation (MFA)'^ support this stand- 
point predicting an increase of Tc as nl^^rii. 

However, at large ric, the Friedel oscillations of car- 
rier spin polarization become significant so that MFA 
becomes inapplicable. More thorough calculations with 
respect to Friedel oscillations (i.e. beyond MFA) corrob- 
orate above statement and show that in DMS undergoing 
FM phase transition Uc cannot exceed some critical value 
related to n^. The reason for that is the oscillations of the 

— 1/3 

RKKY interaction at the scale of l/kp ^ nc , which 
makes impossible long-range FM correlations if rfc^? > 1, 
(r is an inter-ion mean distance)^. 

To properly account for above Friedel oscillations, 
which is indeed a spatial dispersion of inter-ion inter- 
action, we developed so-called random field approxima- 



tion (RFA) in Ref 4. In that work, the ion-ion interac- 
tion has been considered in the context of Ising model. 
This model can be applied for axially symmetric semi- 
conductors with magnetic ions interacting indirectly via 
holes (i.e. with RKKY interaction). The reason for Ising 
model usage for RKKY interaction is a complex structure 
of valence band in DMS that picks out the hole angular 
momenta projections Jz = ±3/2 for lowest heavy hole 
subbands in crystals with distorted cubic or uniaxial lat- 
tice. 

Additionally to spatial dispersion of inter-ion interac- 
tion, there are fluctuations of local magnetic fleld direc- 
tion that also cannot be described in terms of MFA. This 
effect stems from the contribution of transversal spin 
components in effective Heisenberg - like Hamiltonian 
of spin-spin interaction. Namely, RKKY-interaction in 
p-doped undistorted DMS with cubic lattice represents 
this situation. 

In the present paper we pay attention to the fact, that 
transition from Heisenberg spin-spin interaction to Ising 
one, i.e. exclusion of transversal spin components from 
Heisenberg Hamiltonian decreases the system entropy 
and therefore can enhance Tc- For quantitative descrip- 
tion of this effect we present a comparative analysis of 
the RFA-theories for critical temperature Tc in Heisen- 
berg and Ising models. This analysis has been made 
to determine the role of directional fluctuations (inher- 
ent to Heisenberg model) of localized spins in a random 
magnetic field. The uniaxial stresses in typical DMS- 
structures grown on a substrate with some mismatch of 
lattice constants is shown to be the factor responsible 
for appearance of Ising-like interaction between magnetic 
ions spins. Thus, we predict a new effect implying that 
increasing of Tc can be controlled within certain bounds 
by the uniaxial stresses in DMS. In other words, we ex- 
pect that strain engineering can efficiently control the 
value of ferromagnetic phase transition temperature re- 
sulting from the hole-mediated exchange interaction be- 
tween magnetic ions in DMS. 
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II. THEORETICAL BACKGROUND 
A. Heisenberg model 

The Hamiltonian of Heisenberg model for DMS reads 

3<j' 3 

where external magnetic field Hq and interaction J{rjji) 
is measured in energy units (i.e. g^i = 1, ^ is the Bohr 
magneton). The transition to Ising model means, here- 
after, keeping only SjzSj'z - components in scalar prod- 
uct SjSj>. 

The Hamiltonian Q incorporates two kinds of ran- 
domness. First, (the spatial disorder) is that spin can 
be randomly present or absent in the specific j-th cite 
of a host semiconductor. Second, (the thermal disorder) 
is a random quantum state of a spin in j-th cite. These 
spatial and thermal fluctuations can be taken into con- 
sideration by introduction of random field rather than 
mean field. 

In the random field approximation, wc consider every 
spin Sj as a source of fluctuating (random) field 



-J{n-fj)Sj 



(2) 



affecting other spin at the sites fi. In other words, ev- 
ery spin is subjected to some random (rather then mean) 
field, created by all other spins. So, all thermodynamic 
properties of the system will be determined by the distri- 
bution function f{Hr) of the random field Hr. Namely, 
any spin dependent macroscopic quantity (like magneti- 
zation) << A>> reads 



(3) 



<<^>>= / < A>j} f{Hr)dHr, 



where 



<A>..^i^i^^Pll^^ (4) 



Trexp(-7iz/T) 



is a single particle thermal average with temperature T 
and effective Zeeman Hamiltonian Ti.z = H^S. 
The distribution function f{Hr) is defined as 



Latter assumptions with respect to spectral represen- 
tation of 5 function permit to transform Eq. Q to the 
non-linear integral equation for f{Hr) = f{H) in ther- 
modynamic limit. Introducing the probability ni{r)d?r 
for small volume (P'r to be occupied by a particle, we 
obtain 



/(#) - j exp yif{H -Ha)-Q{T) 

g{T)^ [ «'(f)n,(f)dV, 
Jv 

-^(r) =« 1 - exp[iJ(f)^T] >>= 

W{H)fiH)d^H, 
W{H) =< 1 - eyip[iJ{r}ST] >f} . 



(fir 
(2^ 



(6a) 
(6b) 

(6c) 
(6d) 



Eqs (0 represent the integral equation for distribution 
function f{H). In general case this equation can be 
solved only numerically. 

However, in many cases (e.g. for Tc or magnetization 
calculations) it is possible to avoid the solution of the in- 
tegral equation since in these cases it is exactly reducible 
to the set of transcendental equations for macroscopic 
quantities like << 5" >>, 1 < n < 25 of the sys- 
tem. Simplest situation corresponds to the case S = 1/2, 
where only magnetization A4 = — .g/i << S >> (or in di- 
mensionless units m = — 2<<S'>>;(7is g- factor of a 
magnetic ion) is a unique order parameter to be found. 

Thus, in the case of Hq ~ and S — {a are the 
Pauli matrices) Eq. H6a|) takes the form 









1 exp 


ifH 


X 



where 



m= / |tanh^/(iJ)d3if; 



(8) 



To{x) = / n{r) [1 - cos { J{r)x)] dPr; (9) 
Jv 



I{Hr 



5 Hr 



E 



Jin 



Ho 



(5) 



where the bar means averaging over spatial disorder. Our 
RFA approach is based on micro-canonical statistical the- 
ory of magnetic resonance line shape>^ Latter theory as- 
sumes the additivity of local molecular field contributions 
Hr — J2j Hri,j of each particle j (Eq. ((2Jl) as well as the 
non-correlative spatial distributions of magnetic ions. 



Ti{x) — / n(r) sin (J(r')x) d r. 
Jv 



(10) 



Here r = |r| = ^Jt^ + + ^I- To derive (0) we have 
used following property of Pauli matrices: 

exp((T6) = cosh6 -I- ^ sinh6, b=\b\. (11) 
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Substitution of equation lO into Eq.® results in a 
single closed equation for order parameter fa: 



H 



o ; d^H / cPt^ tanh ^ x 
3 / I H 2T 



(2^) 



where 



B = H -m 



(12) 



(13) 



To simplify the vector equation H12|) . we scalarwise 
multiply its left- and right-hand sides by m and inte- 
grate the resulting equation for scalar quantity m = 
over the angle between vectors H and to. The final result 
(see Appendix A for details of its derivation) reads 



TO 



POO 

-6 / 6f/2(t)e- 
Jo 



:Fo{t/2T) 



xi?2 I mTi 



2T 



dt. 



(14) 



S^2(i) = ^csch7ri(l -f ■Ktcothnt), (15) 



Rn{x) 



X COS X — sm X 



(16) 



Trivial solution m = of the equation (|14() corresponds 
to paramagnetic phase. Under certain system parame- 
ters and temperatures, the equation H14fl has nontrivial 
solution that determines the phase transition to the state 
with spontaneous magnetization. 

To find the critical temperature Tc, we use the Landau 
theory with to as an order parameter. For this purpose, 
we may derive (see Appendix B for details of derivation) 
the free energy of the system in the form 



1 f°° 
sin (to.Fi iw)) ^^ 



-ro(t/2T) 



(17) 



In the vicinity of Tc, the free energy H17|l can be substi- 
tuted by a Landau expansion 



Fo 



2A") + —m^A 



H. 

3 ' 



(18) 



2T 



where Fq is a system free energy in a paramagnetic phase. 

It should be noted here that contrary to conventional 
phenomenological Landau expansions of a free energy, 
the coefficients in the function (|18|l have been de- 
rived microscopically within our RFA approach. Free 



energy functions (|17(l and (|18|) give a possibility to de- 
scribe the experimentally observable equilibrium thermo- 
dynamic characteristics (like magnetic susceptibility, spe- 
cific heat etc) of the DMS both in paramagnetic and in 
ferromagnetic phases. 

According to Landau theory of phase transitions, the 
phase transition temperature Tc is reached, when coef- 
ficient 1 - 2Af = in Eq. ((THJ. This is because 
is defined as a temperature, where nonzero infinitesimal 
magnetization appears. Obviously, the same equation 
can be obtained from the Eq. (|14|1 for magnetization in 
the limit to ^ 0. The explicit form of the equation for 
Tr = r5 reads 



-^o(t/2Tf)^, 



(19) 



where = 1/2 in our case. 

Actually, the Eq. H19() determines the T^^ as an im- 
plicit function of system parameters (like rii, Uc etc). 
This function can be considered as a phase diagram that 
separates the region of parameters where the ferromag- 
netic phase with m ^ exists from that where m = 0. 
Latter phase may be paramagnetic or spin glass phase. 
In principle, our RFA method permits to investigate this 
question. This study, however, is beyond the scope of the 
present paper. 

The limit ^ in (P|) gives the relation between 
parameters of the system, which determines the condition 
for ferromagnetic ordering to occur in DMS at T = 0. 
The explicit form of this condition reads 



1 < 



Stt 



^1 (x) 



e-^"^''Ux. 



(20) 



B. Ising model 

Let us consider now Ising model. In this case all effec- 
tive magnetic fields are directed along OZ axis, so that 
scalar product reduces to Sf = Szt, Sz = —1/2,1/2 
and Ea. H6a|) becomes 



f{H) = — / e^"^-^^^^dT, (21) 

Gir) = (^(l- n,(f) (e-'^'^«^-) d'f^'j = 
= .fo(0+iTO.Fi(0, (22) 
where definition of to is similar to Eq.©: 

m = J tanh(^^^ f{H)dH. (23) 

Multiplying Eq.|2U by tanh(^f/2r) and integrating over 
H, we obtain the transcendental equation for order pa- 
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rameter m. The explicit form of this equation reads 



X sin ( mTi { —] ]dt, 



(24a) 



where 



1 f°° / h\ 1 
/ofa;) — — / tanh ( — 1 sin (xh) dh — — — 



TTX 



(24b) 

The equations for free energy and critical temperature 
can be obtained similarly to those for Heisenberg model. 
They read 



F' ^ Fo + _ J B{/^{iit) exp (^-To (^^ 
1 — cos [niTi (2^)) 



-dt. 



F[ = Fo + l^rn^l-A{) + ^m^Ai, 



2T 



2T 



Similar to Eq. (|19|l . the explicit form of equation for 



Tc = T/ reads 



l-2^°°^f/,(t)^,(^)exp(-^„(^))d, 

(25) 

while the equation for FMm region in the phase diagram 
at T = has following form 



°° rl-r 



2 f°° 
1< - / J'i{x)e 
TT Jo 



(26) 



III. DISCUSSION OF RKKY INTERACTION 

A. Non Gaussian fluctuations for 5 = 1/2 

Let us analyze the equations for critical temperatures 
for Heisenberg (|19|l and Ising (|25l) models in more details. 
The difference between them consists only in the form of 
kernels of integrals for Heisenberg (Eq. (|15|l ) and Ising 
(Eq. Ij24b|) 'l cases, so that equation for critical tempera- 
ture in Heisenberg model can be transformed to that in 
Ising model (and vice versa) by replacement of B^i^{t) 

with B[f^{t). 

We start the analysis of these equations from their 
MFA asymptotics. To get this asymptotics, the func- 
tions and in the Eqs ®, ^ should be 
expanded up to linear terms: ^ 0, ^J, 
where J = Jy ni{f)J{r)d^r. After some algebra, the lat- 
ter approximation allows to reduce the Eas H19|) . (|25|) to 



the expressions for critical temperatures in Heisenberg 
T^^'iM = H) or Ising (M = /) models: 

T^^ = 2jj^ B^/,{t)tdt^-J, (27) 

One can see that latter expression is identically the same 
to well-known MFA result for = 1/2 



T, 



MF 



1 



S'(S'+1) / n.,{r)J{r)d^r. 



(28) 



The Eqs ^ demonstrate also that MFA is in- 

dependent of the choice of Heisenberg or Ising model, 

rpMF _ rpMF _ rpMF 

— — -'^c 

Next terms of expansion of the Eqs ©, H1U|I corre- 
spond to Gaussian asymptotics for distribution function 
of local fields. The purpose of subsequent analysis is to 
compare the (actual, i.e. non Gaussian) fluctuations of 
longitudinal components of random field with those of 
transversal ones. 

Since our theory permits to find the distribution func- 
tion f{H) when a spatial dependence of J{f) is assigned, 
we should specify a magnetic interaction in the sys- 
tem. Usually in the problems of carrier-induced ferro- 
magnetism in DMS, the RKKY interaction^ is considered 
as an effective spin-spin exchange interaction resulting 
in FM ordering. To clarify the role of transversal spin 
fluctuations, here we use the simplest possible form of 
the interaction and neglect all possible factors that can 
influence J(r) (such as nonparab olicity of carrier disper- 
sion law etc, see Ref aTllSllQllClllI for more details). Also, 
the stresses may change the form of J(r), see below for 
discussion. 

In the case of simple one band carrier structure, the 
RKKY interaction reads 



J(r) 



-Joxl/^Ri{2kFi 



(29) 



where Xe = ric/^o, ^0 = (|) J^^.^i^/^md, Jc* is a 
carrier-ion exchange constant, TVq = 1/ri is a concen- 
tration of the cation cites, is the density of states 
effective mass. Note that in our single band approxima- 
tion, the effects of stress may influence on Xe and kp^ 
see^ for details. The threshold temperature of ferromag- 
netic ordering in MFA now can be found by evaluation 
of integrals H28|l with respect to (|29|l : 



4/3 



4/3^1/3^ 



(30) 



Here, the factor Jqx^ is independent of carrier concen- 
tration, Xi — riiQ is a molar fraction of the magnetic ions. 
The ratio of electron and magnetic ion concentrations 
v = ric/ni = Xe/xi plays a crucial role in our theory as 
a parameter separating the cases of relatively small fluc- 
tuations with V <^ Uc and that of large ones with v ^ Vc', 
parameter Vf. is indeed a dimensionless critical concentra- 
tion (corresponding to equality sign in expressions H2U|) 
and H2t)|l for Heisenberg and Ising models respectively). 
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The functions ^o(a^) and Ti{x) in Eqs.lO, ifTUIl (with 
respect to Ea. (|29|l ') assume foUowing form in the case of 
homogeneous magnetic ions distribution, ni{r) — rii — 
const, 

foiO = / {I- cos {^R4y))}y'dy, (31) 



(32) 



(^i(0 = / sin {^R4y))y^dy. 
Jo 



In the case of spin S = 1/2, the resuh ^ of MFA 
can be recovered from Eqs H31|l . I|32|l by their expan- 
sion up to hnear terms tpo 0; 9?i (0 f (of 
course, with their further substitution into Eq. I|14|l or 
Eq. H24a|l '). Gaussian asymptotics of distribution func- 
tion corresponds to the next term of expansion of the Eq. 

To account for real (non-Gaussian) distribution of fluc- 
tuating local field, we do not expand Eqs and 
calculate them numerically. 

In dimensionless variables, the equations for critical 
temperatures for both above models assume following 
form 



1 

^1/2 



— oo 
t 
hi 



t 



= exp 



26m 
1 



^1/2 
/ t 



dt = 1, 



(33) 



where 9m = T^' / (Joxt^^)', M stands for H (Heisenberg 
model) or / (Ising model). The result of calculation of 
T^^/T^^^ with the help of Eq. ijSSI) as a function of v is 
reported in the Fig. la. It is seen, that there are curves 
that separate the areas of system parameters (including 
temperature) where FM or non-FM phases occur. 

Our results show the limited area of concentrations 
Q < V < Vc which allow FM ordering in both considered 
models. So, we have found — 0.2473 for the Ising 
model and 0.0989 for the Heisenberg model. Our results 
also show the critical character of dependence Tc = Tc{v) 
that can be well approximated by the function 



Tr ~ T, 



MF 



(1 - y/^cY 



2Att 



[l-y/ycT (34) 



with A ~ 0.47 and A ~ 0.63 for Heisenberg and Ising 
models respectively. 

Thus, the fluctuations of transversal spin components 
suppress a tendency towards FM ordering in the range 
of concentration ratios < < 0.099. Moreover, in the 
range 0.099 < v < 0.247, our results predict impossibility 
of FMm in DMS with Heiscnberg-like spin-spin interac- 
tion while in this interval of v FMm can still occur in 
DMS with Ising- like interaction between spins. Since 
Ising model is inherent to uniaxially stressed semicon- 
ductors, latter conclusion means that uniaxial distortion 
can effectively inhibit transversal spin fluctuations thus 
enhancing Tc- 
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FIG. 1: The phase diagram of the systems under considera- 
tion, a) general (non-Gaussian) case for spin 1/2, b) Gaussian 
approximation for spin 5/2. H - Heisenberg model, I - Ising 
model. Region 1 corresponds to FM state for both models, 
region 2 corresponds to FM state for Ising model and non FM 
state for Heisenberg model, region 3 corresponds to non FM 
state for both models. 



B. Gaussian fluctuations for S — 5/2 

The equation (|19|l for critical temperature Tc is exact 
in the framework of our model, which means that it incor- 
porates non-Gaussian fluctuations. At the same time it 
is restricted by the case of ion spins S = 1/2. Practically 
important case of Mn?^ ions with 5 = 5/2 needs a spe- 
cial consideration. Mathematically, the problems arise 
for S > 1/2, when exponential function of spin operator 
is no more a linear function so that Eq. is no more 
valid. In this case we should use Sylvester theorem for ex- 
ponent of aHermitian matrix (25-1-1) x (2S'-|-l)(see, e.gi^^ 
for more details). Appendix C presents the final result 
of such calculations for W{H) and S = 5/2 (Eq.j6dJ). 
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It is apparent from Appendix C that the case S = 5/2 
involves much more order parameters (such as a^.n = 
{c'^Vn) = / d^iinf{H)(fH, m = 0,2,4; n = 0,1,2, and 
hm,n = (c"'0„), m = 1,3,5; n = 0,1,2) than in the 
case S = 1/2. Not all of them play a crucial role in 
a formation of FM phase in DMS. For example, if the 
random field fluctuations are almost Gaussian, we can 
find the expectation values of spin operators Q with the 
aid of only two parameters, representing first and second 



moments of the distribution function f{H). 

Our calculations show that Gaussian approximation 
for f{H) adequately describes the actual phase diagram 
of DMS (except for close vicinity of critical concentration 
which is not practically important). That is why the 
analysis of the case of 5 = 5/2 can be performed with 
sufficient accuracy in Gaussian approximation. To do 
this, the Eq. (|G6|I should be substituted by its expansion 
int = J{r)T up to the second order. The result reads 



J 



W{H) = ic 



n , 3iJ 1 , H 

3 coth — - 2 coth — 



J{f)i 



l + 4c^ - 



(3c2 - 1) (l + 4coshf ) 
1 + 2 cosh ^ 



(J(f)r)^ 



(35) 



r 



We apply this result to obtain the equation for critical 
temperature. The magnetization and mean value of ran- 
dom field H are negligibly small at this temperature that 
suggests that f{H) is isotropic. Thus, performing angu- 
lar averaging over the directions of H in Eq. H6c|l , we can 
put = 1/3 in the Eq. (jSHJ- 

To relate this result to the case S = 1/2, we note that 
first term of Eq. H35|l includes Brillouin function for spin 
S = 5/2 

„ fSH\ , 1 , , (5 + l/2)iJ 1 , B 
Bs\ = IH ) coth -^^ coth — . 

(36) 

After simple algebra, the Eq. takes following form for 
arbitrary spin S 



= i-msj{f) + l-sis + 

TO 2 



where 



H fSH 



fiH)d^H. 



(37) 



(38) 



One can see that the above equations are formally sim- 
ilar to those for S* = 1/2 if we expand Eqs © and HlOfl 
up to first nonvanishing terms. The Eq. (|37|l determines 
the components of Fourier image of distribution function 
in Gaussian approximation 



»(^) 



n(f)J"(f)dV. 



5Jt, 



(39) 



Equations H38|l and H39|) reduce the problem of de- 
termination in Gaussian approximation for arbitrary ion 
spin to the case of non-Gaussian fluctuations for spin 
5=1/2 considered above. Namely, after substitution 
of these equations to Eq. H19|) along with proper gen- 
eralization (for the case of arbitrary spin) of the func- 
tion yB^jlOi can use this equation to find Tc for any 
S > 1/2. Calculations performed in the same manner as 
Eqs. and (|X8|l yield 



coth ■ 



TTt 



2S+1 



cothTrt H 




csch' 



2 TTt 

2S+1 



25*+ 1 



csch irt 



(40) 



r 



It can be readily shown that in the case = 1/2 the Eq. 
(|40|l reduces to simpler form (|15|l . 

Substitution of (|39|l to H20() permits to obtain the nec- 
essary condition for FMm formation at zeroth tempera- 
ture in the form (see also Eqs. (pn|l and 



J> 



37r5-|- 1 



S 



J2. 



(41) 



This inequality has clear physical meaning: small disper- 
sion and positive (FM) inter-ion spin-spin interaction fa- 
vors creation of FMm in DMS. Also, large spin is also 
preferable for FMm formation due to relatively small 
(quantum) fluctuations of its transversal components. 

For comparison, we also consider Ising model for 
S > 1/2 in Gaussian approximation (despite the possi- 
bility to account for non-Gaussian fluctuations of ran- 
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dom field for arbitrary S > 1/2).'' Expansion of the 
function G{t) fEg. 1)22(1 ') with arbitrary S and m = 
I^oo (^) f{H)dH up to second order leads to the 
distribution function in the form 121() with G{t) in the 
form ((22(1 . where 

^o(0^M.J^^; ^^Q-SJr. (42) 
Parameters 

IZ^'i^) ^^^^'^^ ^^^^ 

and 

- {{Si)) = 1^ Qs (f ) f{H)dH, (44) 

Q5 (f) -^(^+l)- 
-5coth(|,)B5(^), (45) 

should be found self-consistently with respect to Eqs. 
((21(1 . ( 122(1 and 1(42(1 . Thus, in Gaussian approxima- 
tion, Ising model needs self-consistent determination of 
two parameters, {{Sz)) and {{S^)), while in the case of 
Heisenberg model single parameter << \S\ >> needs to 
be determined. This is because for Heisenberg-like in- 
teraction << 5^ >>= S{S -I- 1) = const (Eq. (jSHJ) so 
that this parameter depend neither on temperature nor 
on field distribution. 

The critical temperature can be found from the Eqs 
(gni and ijlll) as TO — > 



V7rM2J2 

"io '''''' [WTTJl^)' ''-^ ^''^ 

where 9 = T^/T^^^ is a ratio of actual critical tem- 
perature and that obtained in MFA approximation (Eq. 
((28(l ). In this approximation, the necessary condition to 
form FM state at T = in DMS has the form of following 
inequality 



J>^-J2. (48) 

The independence of Eq. 1(48(1 of spin can be thought of 
(see Eq. ((41(1 for comparison) as a lack of transversal spin 
components contribution to random field in Ising model. 

Now we are ready to compare Heisenberg and Ising 
models for specific case of RKKY interaction ((29(1 . In 



Gaussian approximation we have to evaluate only two 
integrals 

7^^Jo4/v/^ J^=^.(7)^ (49) 

OTT 5 

Note that the equations 1(49(1 demonstrate the relation 
between Friedel oscillations and first (cx J) and second 
((X J^) moments of distribution function of random mag- 
netic fields. Namely, at frequent Friedel oscillations (i.e. 
large ric and v) the influence of dispersion [J"^)^/"^ (which 
is "responsible" for disorder in the system) prevail over 
trend to order the system due to mean value Jj^ This 
means that frequent Friedel oscillations at the scale of 
mean inter-ion distance inhibit ferromagnetism in DMS. 

Substitution of gSjl into gSJ , (gTJ permits to find the 
dependence 0{v) for practically important case S — 5/2 
(see Fig. lb). One can see that qualitatively situations 
for S = 5/2 and 5=1/2 are similar: the region of FM 
state is significantly larger for Ising model than that for 
Heisenberg model. But there is also a difference. Com- 
parison of Fig. la and lb shows that the area between 
curves Tc{nc/ni) for Heisenberg (H) and Ising (I) models 
at low temperatures is smaller for S — 5/2. This means 
that for S — 5/2 quantum fluctuations (which is the only 
possible fluctuations at T = 0) are inhibited as compared 
to the case S — 1/2. 

This is actually a reflection of the well known fact that 
the larger the magnitude of the quantum number of the 
spin, the "more classical" it is, i.e. the smaller is the 
contribution of quantum fluctuations of its transversal 
components. At T ^ additional thermal fluctuations 
appear. It is seen from the Figure, that at T < the 
extension of FM phase due to enhancement of Tc is al- 
most the same for S — 5/2 and S = 1/2. This means 
that the thermal fluctuations of the spin do not sensitive 
to its value. 

IV. CONCLUSIONS 

In this paper we have proposed a new mechanism of the 
enhancement of FM phase transition temperature Tc by 
the uniaxial distortion of DMS. This prediction is based 
on comparative analysis of Heisenberg (inherent to undis- 
torted DMS) and Ising (inherent to uniaxially distorted 
DMS) models. The analysis of above models has been 
carried out in the framework of our recently developed 
formalism^ random field method. This method, which 
can be regarded as a substitution of conventional MFA 
for disordered systems with given J(r), permits to derive 
self-consistently the equations for order parameter to and 
the free energy functions of DMS. 

Now we discuss in more details the infiuence of stresses 
on magnetic spin Hamiltonian of DMS 1^. We consider 
such influence in two steps. In first step we consider the 
influence of the stress on the effective spin-spin interac- 
tion potential and in second step we consider the operator 
part of the Hamiltonian. 
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Since the effective potential of indirect spin-spin inter- 
action strongly depends on the band structure of specific 
semiconductor sample, which is affected by the stress, 
this potential by itself may also depend on stress. This 
influence^ is manifested via both density of states (in our 
case we use the effective mass of density of states nid (I29|l ) 
at the Fermi level and concentration of free holes, related 
to pinning of the Fermi level by defects and impurities of 
different nature. As it was noted in Ref.^, the influence of 
pressure on the density of states is small. The influence 
via pinning centers depends entirely on their nature. We 
can imagine the situation when the influence of pressure 
on the concentration is also small. For example, in the 
cases of the absence (or small number) of pinning cen- 
ters or " synchronous movement" of pinning centers with 
the valence band edge shift due to deformation potential, 
this effect is negligible and our mechanism of influence 
of the pressure will be decisive. Here we would like to 
emphasize that there are no general problems to incor- 
porate the possible dependence of the concentration on 
the stress into our consideration (this is simply one more 
modification of J{rjji) in (^). If we do so, the consid- 
ered effect of elimination of transversal spin components 
by the stress, which has not been discussed in the litera- 
ture, is an additional factor enhancing Tc. 

To discuss the stress effect on the operator part, we 
note that in the typical case of p-type DMS with cubic 
lattice the transversal spin-spin interaction can be sub- 
stantially inhibited by the uniaxial stress of a crystal. 
Really, such stress splits the valence band edge to light 
and heavy hole subbands. The resulting heavy hole (HH) 
states are characterized by angular momentum projec- 
tion Jz = ±3/2. For such states, the spin-flip scattering 
of these holes on magnetic ions is forbidden due to an- 
gular momentum conservation. Hence only longitudinal 
(along the distortion axis) Z-components of the spins will 
be present in the resulting indirect spin-spin interaction 
via above heavy holes. 

However, aforementioned spin-flip processes are not 
forbidden both for light hole (LH) states (with spin 
projection ±1/2) and for transitions between LH and 
HH states thus contributing to the indirect inter- 
action of transversal spin components of magnetic 
ions. Hence, the deformational splitting of a va- 
lence band edge, leading to preferential occupancy 
of the heavy hole subbands gives the anisotropy 
of indirect spin-spin interaction in the form Ti = 
Si>i [■J\\i,j^ZiSzj + J±ij {SxiSxj + SyiSyj)] ■ The ra- 
tio 7 = J±i,j/J\\i,j (0 < 7 < 1) should be monotonically 
decreasing function of the ratio of heavy holes concentra- 
tion hhh to their complete concentration nn- Thus, if 
deformational HH - LH splitting exceeds the Fermi en- 

I 



(2Trf J J H 



ergy of the holes, the ratio nHH/nn reaches its maximal 
value {nHn/nH 1), corresponding to Ising Hamilto- 
nian (7 — )■ 0) of the indirect spin-spin interaction. In 
this case the region of ferromagnetic state of DMS ex- 
pands substantially towards both higher carriers (holes) 
concentration and higher temperatures. 

The detailed theoretical description of all above ef- 
fects, which is intimately related to the parameters of 
specific DMS sample and experimental conditions, can 
be developed within the framework of presented theory 
for any particular case. However, such calculations can 
be done only numerically. Note, that the magnitude of 
effect which we predicted is very sensitive to the holes 
concentration Uh. The problem of its correct determi- 
nation, to the best of our knowledge, is related to the 
influence of anomalous Hall effect. Latter effect lowers 
reliability of Uh determination from Hall effect datai^^. 

For better illustration of our effect, we estimate now 
the magnitude of Tc increase for typical ferromagnetic 
DMS Gai-xMrixAs deposited on GaAs and GaP sub- 
strates. The mismatch Aa of lattice constant a leads 
to biaxial strain that splits the valence band with de- 
formation potential b — —1.7 eV by the value 5Ei,2 — 
2|6e^^|ji^ where zz- component of strain tensor e^^ = 
— 2(Aa/a)ci2/cii, and the ratio of elastic moduli in 
GaAs C12/C11 = 0.453. For x — 0.035 (or concentra- 
tion rii — 7.76 • 10^° cm""^) , the relative mismatch 
|Aa/a| = 0.002 for GaAs substrate and |Aa/a| = 0.036 
for GaP substrate.^ We can see that for GaAs substrate 
small valence band splitting 5Ei^2 — 6 meV cannot sup- 
press interaction of transversal Mn-spin components for 
typical concentration ric = 10^'^ cm~'^ which corresponds 
to Fermi energy ep ~ 80 meV, whereas for GaP substrate 
5Ei^2 — 109 meV > sp- Thus, our mechanism predicts 
the enhancement of Tc for Gai-xMn^As on GaP sub- 
strate by the factor 1.64 (see Fig. lb for v — nc/rii — 0.13) 
as compared to the same DMS but on GaAs substrate. 

Let us finally note that in the present paper we con- 
sidered the enhancement of Tc due to RKKY interaction 
only. But there are also other mechanisms, which can 
lead to appearance of ferromagnetism in DMS, see^ for 
details. These mechanisms will be eventually reduced 
to the Hamiltonian with modified potential J{fj,j>). 
Thus for quantitative discussion of these mechanisms it 
is sufficient to substitute the corresponding modified po- 
tential to our self-consistent equations. 



APPENDIX A 

We begin with equation ()12|) for magnetization. 



— exp{-.Fo(^) +2fs}, (Al) 
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First, we pass to spherical system for f and integrate the Ea. ljAl|l over f directions. This yields 

sin Bt 



1 f H H 

^Jd^Hl r^rfr-tanh-exp 



-^0 



Bt 



B = \B\ = \J m -2cos6mHTi{T/2)/T + {mTi{T/2)/T)\ 

where 9 is an angle between vectors rh and H . 

Next step is a scalar multiplication of its both sides by fh that yields 



(A2) 
(A3) 



= ^-^ / d^H j r'^dr^^^ tanh ^ cxp 



Hrh 
IT 



H 
2T 



or 



m= — - / (PH / r^dr cos tanh — exp 
27r^ ./ In 2T 



sin Bt 
Bt 



sin Bt 



Bt 



(A4) 



The rotational invariance of scalar product permits to point B along z axis and integrate over angular variables in 
[4|). This yields 



1 

2^ 



dip I sin 9d9 I H'^dnl r^dr cos 6* tanh ^ exp [-J^o (^) 

.In .In .In 21 L \ Z / 



sin Bt 



where B is defined by HA3|I . Consider 



Change of variables 



1 



27r 



/= — / dip f sin.9 cos 9^^^^-^^ d9. 

T ./o .In B 







cos V = z 



+ {mJ^^{T/2)/Tf ~B^ 
2mHTi{Tl2)lT 



with the help of IjASp reduces it to the form 



/ = 



T {mH Ti{t) / t) J H+mJ^i[T)/T 
= ATlR2{mTi{T/2))R2{HT), Rn{x) 

Substitution of IjASfl into (|A4|1 gives 



H-' + {mTi{T/2)/Tf -B^ 



sin BTdB 



X cos X — sm X 



2 /"OO pCG 

m — — dH I dT tanh — exp 
7r7o Jo 2r 



-^0 



i?2 imHiiT/2))RoiHT). 



(A5) 



(A6) 



It is also possible to integrate over H in (|A6|) 

/ tanh — Ro{Ht) dH = -dnTB^Ut), t = rT. 

^0 2r 



Sf/2W = 



sinh 7rt + irt cosh 7ri 
3 sinh^ Tit 



(A7) 



(A8) 



APPENDIX B 

We start the derivation of the free energy from the 
equation p4|) for magnetization (order parameter). We 
rewrite it in the form 



t 



m 



6^°°6f/2(^i)e-^"^*/'^^i?2 (mT^ ] ] dt ^ 0. 



(Bl) 

With respect to substitution tT = t this gives final equa- Now we recollect that if we have a free energy F of a 



tion for m in the form (|14|) of the text. 



system, then the equation for order parameter (in our 
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case Eq. (HU) should minimize it. In other words, Eq. Condition ljB2|) is a simple differential equation for F, its 
(|ijlll should be equivalent to condition solution yields 

dF 



dt 



This is indeed equation H17() from the text. 

To get Landau expansion of ljB3p . we simply expand sin (rn!Fi {■^)) in Taylor series at small m. This gives 



F ^ Fo + Gj^ dt + - "^'^f + ^^"'^3'' + (B4) 



Paying attention that second term in l|B4(l does not de- 
pend on m and hence just renormalizes Fq^ we easily 
obtain Eq. 118|) from the text. 



APPENDIX C 

We are looking for the expression Ij6d|l for the case 
S — 5/2. We introduce notation t — J{r)f and assume 
T = 1 (i.e. H means H/T). In these notations Ea.ll6c 
reads 



We introduce the functions 



W{H) = 1 



Tr |exp 


iSt 


exp [-SH^ 1 




exp 





(CI) 



The denominator in Ea. ljCl|l can be immediately evalu- 
ated in a reference frame rotating around quantization 
axis 



Tr|exp -SH | 



sinh 3H 
sinhiJ/2' 



(C2) 



cosh (nH) 
" 1 + 2 cosh {2H) ' 

cosh {nH) (coshiJ — 1) 
^" " sinhSi? ^ 

cosh (niJ) sinhi7/2 
cosh iJ/2 + cosh 3H/2 + cosh 5H/2 



(C3) 
(C4) 



and cosine of the angle between vectors H and f \\ t, 



cos { H,f] = 



Ht 



Ht 



(C5) 



In these notations after lengthy calculations we can get 
the expression for the trace in numerator of Eq. IjCljl . 
With respect to Eq. (|C2|) . this expression assumes the 
form 
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W{H) = 1 - ^cos^ |(3- 14c2 + 15c'') (/)o -4(1 -6c2 + 5c4) 01 +5 (1 - c2)^02} - 
1 3t 

-- cos — { (-1 + 38c2 - 45c'') </>o + 4 (1 - 12c2 + 15c^) (t>i + 5 (l + 2c^ - Sc") 02 } - 
8 2 

1 5t 

-- cos - { (3 - lOc^ + 15c*) 00 + 4 (l - 5c"') 0i + (l + lOc^ + 5c*) 02} + 

o Z 

+ J sin ^ { (3 - 14c2 + 15c*) 7?o - 4 (l - 6c^ + 5c*) r/i + 5 (l - c2)%2} - 
ic 3t 

-- sin — {(3 - 26c^ + 15c*) r/o - 4 (3 - 4c^ + 5c*) Tyi - 5 (3 - - c*) 772} + 
8 2 

+ ^ sin ^ {(15 - lOc^ + 3c*) r,o + 4 (5 - c*) ryi + (5 + lOc^ + c*) 772} . (C6) 
8 2 
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